Topologically induced avoided band crossing in an optical chequerboard lattice 
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We report on the condensation of bosons in the 4th band of an optical chequerboard lattice 
providing a topologically induced avoided band crossing involving the 2nd, 3rd, and 4th Bloch bands. 
When the condensate is slowly tuned through the avoided crossing, accelerated band relaxation 
arises and the zero momentum approximately C4-invariant condensate wave function acquires finite 
momentum order and reduced C2 symmetry. For faster tuning Landau- Zener oscillations between 
different superfiuid orders arise, which are used to characterize the avoided crossing. 

PACS numbers: 03.75.Lm, 03.75.Hh, 03.75.Nt 



Topologically non-trivial band structures provide the 
basis of intriguing forms of quantum matter such as high- 
temperature superconductors, unconventional supercon- 
ductors in heavy-fermion compounds, quantum Hall sys- 
tems or the recently discovered topological semi-metals 
in semiconductors with strong spin-orbit coupling [l]-iij. 
Topological classification of band structures and band 
crossing points has been a topic of vibrant theoretical 
research recently The unique possibilities to p re- 

pare tailored periodic potentials in optical lattices 11 , 12j 
have recently inspired numerous proposals t o use them 
for clean simulations of topological matter 
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the lowest band artificial gauge fields may permit to engi- 
neer non-trivial topological features, however, at the cost 
of significant experimental complexity fljj . In the higher 
bands of two- or three-dimensional non-separable band 
structures topological properties can be readily obtained 
by an appropriate choice of the lattice symmetries. The- 
oretical work has predicted that band relaxation can be 
kep t moderate for bosonic atoms excited to higher bands 
MEI and recently condensates in the P and .F-bands 



of a square optical lattice were reported [22|, [23j and in- 
vestigated theoretically [24| . A complementary effort in 
excitorr-polariton systems has led to a D-band conden- 
sate 



25|. 



In this work, we combine the tunability of optical lat- 
tice potentials with targeted condensate formation in an 
excited band to implement and investigate a topologically 
induced avoided band crossing. Extending techniques de- 



scribed previously [22j, |23j , we produce an optical lattice 
potential with two classes of sites (denoted as A and B- 
sites in Fig. 1(a)) arranged as the black and white fields 
of a chequerboard. In the xy-plane the optical potential 
with tunable parameters Vb > and AV is given by 



V(x,y,e) = -^|r? (e ik * 
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FIG. 1: a, The lattice comprises two classes of sites with 
different well depth denoted by A and B. The grey area de- 
notes the Wigner Seitz unit cell, b, The first BZ with a path 
marked connecting the F-, X- and M-points, respectively, c, 
Energies of the 2nd (B 2 ), 3rd (B 3 ), and 4th (B 4 ) Bloch bands 
plotted versus AV for a lattice with perfect C4-symmetry for 
Vb = 7.0-Ercc and AV C = 6.lE 1C c. The detail on the lower 
edge shows a plot of these bands for AV" = AV C within the 
1st BZ along the trajectory illustrated in (b). d, The bands of 
c are replotted for a lattice with weakly broken C4-symmetry 
for Vb = 7.8E IOC and AV C = 6.1 E ICC . The details on the 
lower edge show plots analogue to that shown in c, however, 
for AV = AV C ± E Tec . The (thin, grey) dashed dotted lines 
mark the energy of the 4th band at the F-point 



with e* = (77, e x , e v ) and carefully measured quantities 77 = 
0.98, e x = 0.93, e y = 0.87 (accounting for imperfect optics 
in the experiment). The angle 9 can be adjusted to better 
than 7r/200 by means of an interferometric optical set- 
up [1^]. The effective well depth difference is defined as 
AV(6) ee -V r](l + e x )(l + e y )cos(6), i.e., AV(tt/2) = 0. 
In the z-direction a weak harmonic potential is provided 
giving rise to elongated tubular lattice sites. 

We denote the Bloch bands as B n ,n G {1,2,...} or- 
dered according to increasing energies. In Fig. 1(c) the 
bands B 2 , B 3 and B4 at the T-point (cf. Fig. 1(b)), de- 
rived from a two-dimensional (2D) band structure cal- 
culation for a lattice with Vb = 7.0E ICC and perfect C4 
symmetry (i.e., setting 77 = e x = e y — 1 in Eq.(l)), are 
plotted against AV. Here, E lec = (Hk) 2 /2m = recoil en- 
ergy, m = atomic mass, k = 2ir/\, and A = 1064 nm. 
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The plot shows that at the critical value AV = AV C 
(= 6.1 E ICC for Vq = 7.0 E iec ) all three bands become de- 
generate and the B2- and the £?4-band form a Dirac cone 
in the centre of the first Brillouin zone (BZ) with the 
locally flat £?3-band intersecting its origin. This can be 
seen in the detail on the lower edge of Fig. 1(c) showing 
the involved bands within the first BZ along a trajec- 
tory connecting the points T, X, M, V of Fig. 1(b). For 
AV > AV C (AV < AV C ) the B 4 -band (B 2 -band) sep- 
arates, whereas B2 and B% (B3 and B4) remain degen- 
erate, thus forming a topologically protected quadratic 
band crossing point (TQB) [6j at the centre of the first 
BZ. This structure is robust against changes of AV and 
Vq as long as the CA symmetry is sustained. 

In our experimental realization, due to technical imper- 
fections accounted for by the value of e specified below 
Eq.(l), the C4 symmetry of the lattice is weakly broken. 
The lattice potential may be decomposed into a main 
contribution with perfectly C4-invariant unit cell and 
a small perturbation with no rotation symmetry, which 
acts to lift the degeneracies in Fig. 1(c), thus leading to an 
avoided band crossing with an energy gap Ai?2,4 on the 
order of a small fraction of E rcc , as shown in Fig. 1(d). A 
band calculation shows that 8AE2 4,/dVo < 0.02, i.e., the 
gap AE2.4 ~ 0.26 EVcc is nearly independent of Vq. On 
the left side of the anti-crossing (AV < AV C ) the bands 
B3 and B4 rapidly approach an approximately constant 
separation AE^ « 0.13 -E rcc , which is nearly indepen- 
dent of V {dAE 3A /dV < 0.002). Within the 1st BZ 
the energy AE 3 ^ appears as the gap introduced into the 
TQB at the T-point. An analogue Vo-independent gap 
AE2,3 arises for the bands B2 and B$ on the right side 
of the anti-crossing (Ay > AV C ). The energy gaps have 
been calculated according to the definition AE n rn (Vo) = 
Mm AVe[AVc _2 EtecA v c +2E tec ][E m (AV,Vo) - E n (AV,V )] 
with E n (AV, Vq) denoting the energy of the nth band. 
The robustness of these gaps (with sizes on the order of 
ten nanokelvin) against changes of Vq and AV indicate 
that, despite the broken C-j-symmetry, some topological 
character of the bands appears to be preserved. A theo- 
retical method (extending the Berry flux concept @,[i3l) 
to quantify this residual robustness is yet to be conceived. 
The details on the lower edge of Fig. 1(d) show TXM- 
trajectory plots of the bands within the first BZ away 
from the avoided crossing at AV = AV C ± E rcc . 

As the initial step in our experimental protocol we 
produce a sample of about 10 5 rubidium atoms ( 87 Rb) 
in the F = 2,mp — 2 electronic ground state con- 
densed in the global minimum of the £?4-band for a value 
= AV = 9.0E rcc = AV C + 2.9E 1CC located on the 
right hand side of the avoided crossing in Fig. 1(d), i.e., 
where the £?4-band is well separated from all other bands. 
This value is chosen in order to maximize the decoherence 
time of the condensate fraction to about 26 ms, which 
corresponds to a lifetime of the i?4-band poplation of 45 
ms. For condensate formation (following Refs. [22, 23 1) a 
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FIG. 2: Momentum spectra of the i?4-band condensate are 
shown in a and b for (AV — AV C ) / E lec — 2.9 and in c and d 
for (AV — AV C ) / E ICC = 9.9, respectively. The experimental 
results in a and c are compared to calculations in b and d. 
The AV dependence of the populations of the Bragg peaks 
identified in the dashed black rectangle in (d) is shown in (e) 
(experimental observations) and (f) (calculations). 



deep ground state lattice (Vo ~ 15Ej- cc ) is prepared with 
AV <C such that only the deeper S-sites are occupied 
and tunneling is suppressed. Then, AV is switched to 
AVi, rendering 6-sites more shallow than .A-sites, and 
Vq is reduced to 7.8 E vcc in order to allow for tunneling. 
Within 10 ms the system condenses into the minimum of 
the i?4-band and we may subsequently tune AV to any 
final value of interest. 

Fig. 2 shows that in a wide range of values AV we ob- 
serve momentum spectra with sharp Bragg maxima, re- 
sulting from the coherent condensate fraction residing at 
the T-point of the i?4-band. The panels in (a) and (c) 
show the experimental observations for AV = AVi = 
AV C + 2.9E rcc and AV = AV C + 9.9E IOC , respectively. 
The panels below in (b) and (d) present corresponding 
calculations based upon single particle Bloch-functions 
using the potential of Eq.(l), showing good agreement 
with respect to the relative sizes of the higher order Bragg 
peaks. For AV > AV C , the Bloch- function corresponding 
to the r-point of the i?4-band comprises local ls-orbits 
in the shallow wells and local 3s-orbits (with one radial 
node) in the deep wells. For values AV in the vicinity 
of AV, (as in (a) and (b)) the ls-orbits hold most of the 
population, while the 3s-orbits are only marginally pop- 
ulated. Hence, the envelope of the momentum spectrum 
is mainly due to the Fourier transform of the compara- 
tively dclocalizcd Wannier- function of the ls-orbit, such 
that only a few higher order Bragg peaks appear. Larger 
values of AV yield increased population of the 3s-orbits. 
Because of their larger momentum components, higher 
order Bragg peaks become increasingly visible. This is 
more systematically studied in panel (e) , which shows the 
AV-dependence of the populations of the Bragg peaks 
on the descending diagonal in (d) indicated by the grey 
dashed box. The corresponding calculation in (f) shows 
striking agreement with the observations. For values of 
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FIG. 3: Band mapping plots recorded for (AV — AV C ) / E leG = 
4.5 and holding times 20 [is, 20 ms, 100 ms. In the lower right 
corner a map of the theoretical BZs 1,2,3,4 is plotted. The 
black circles indicate the observed condensation points at the 
energy minima of the i?4-band, which connect to the centre 
of the 1st zone via reciprocal lattice vectors (white arrow), b, 
The decay time of the 4th (blue disks) and the refilling time 
of the 1st band (black diamonds) are plotted versus AV. The 
black line plots show the scaled overlap integrals Ii (solid), 
I2 (dashed), ^3 (dashed dotted) defined in Eq.(2). 

AV > AV C the i?4-band is C4-invariant to very good ap- 
proximation with its global minimum at the T-point. In 
this case the zero momentum Bragg peak is by far the 
most populated. If we tune AV to values below AV C , 
thus passing the avoided band crossing, significant oscil- 
lations of the Bragg peak populations arise, which are 
discussed below. 

In Fig. 3 we study the collision induced relaxation of 
the i?4-band population. After condensate formation, 
AV is adjusted to some desired value within 200/zs and 
the atoms are held in the lattice for a variable time. Fi- 
nally, the lattice potential is adiabatically turned off, and 
the atoms are allowed to expand for 30 ms. This provides 
us with images of momentum space, in which the pop- 
ulation of the nth band is mapped into the nth BZ. In 
(a) an example is shown for (AV — AV c )/E rcc = 4.5, 
where the effect of the avoided crossing is not relevant. 
Initially (for 20 [is holding time), mainly the 4th BZ is 
populated, with a significant fraction of the atoms re- 
siding at the condensation points (±hk,±hk). As seen 
in the pictures for larger holding times, the population 
of the 4th band directly decays into the 1st band (with 
approximately exponential time dependence), while the 
2nd and 3rd BZs are not markedly involved. Only in the 
vicinity of the avoided crossing, a more complex decay 
dynamics arises with the 2nd and 3rd bands initially ac- 
cumulating significant populations before finally the 1st 
BZ is refilled. This becomes visible in (b) (showing the 
relaxation times for the 4th and the 1st bands versus 
AV), where close to the anti-crossing (AV ~ AV C ) the 
decay of the 4th band (blue disks) is faster than the re- 
filling of the 1st band (black diamonds) . The plot shows 
a pronounced resonance around AV = AV; with notably 
long lifetimes above 40 ms. This may be qualitatively 
explained by the observation that around AVj most of 
the atomic population resides in the local ls-orbits of 
the shallow wells, where it is protected from collisional 
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FIG. 4: a, Momentum spectra (left = observation, right 
= calculation) on the left side of the anti-crossing (AV = 
AV C — 1.315 rC c). A sufficiently large gap (~ 0.8E rcc ) was 
adjusted, in order to prevent Landau-Zener dynamics, b, 
Landau-Zener dynamics for a gap of ~ 0.26-E rC c. The vis- 
ibility of the ±(1, — l)ftfc-Bragg peaks (blue disks and green 
triangles) and that of the (0, 0)ftfc-peak (red squares) are plot- 
ted versus the holding time after rapid tuning over the anti- 
crossing to (AV — AV c )/£'rcc = —3.0. c, Fourier spectrum of 
the oscillating (—1, l)fifc-peak with first and second harmonic 
components indicated by black arrows, d, Blue disks denote 
the first harmonic frequencies derived from plots as in (b). 
The solid line shows the calculated energy difference between 
the 4th and the 2nd band. 

decay, because locally no state with lower energy is avail- 
able. This assertion is supported by plotting the integrals 

UAV) = (2) 

\jh d2r Pl Io d2r Pl 

where o denotes the unit cell in configuration space, 
Pn = I0n| 2 an d 4> n ,n G {1,2,3,4} denotes the Bloch- 
function of the nth band for zero quasi-momcntum calcu- 
lated for the optical potential in Eq.(l). These integrals 
measure the spatial correlations of the particle densities 
in the 4th and the nth bands at the T-point. As shown 
in (b) the maximal lifetime of the 4th band arises, where 
these overlap integrals are small. The observation in (a) 
that, apart from the vicinity of the avoided crossing, di- 
rect decay to the 1st band appears as the dominant pro- 
cess, corresponds to the fact that I\ displays the largest 
decrease near the relaxation time maximum. 

As already indicated in Fig. 2, when we drive the con- 
densate across the avoided crossing at AV — AV C , we 
observe momentum spectra with oscillating populations 
of the Bragg maxima. For small band gaps of less than 
about half a recoil energy, the rapid band relaxation oc- 
curring near the avoided crossing (cf. Fig. 3(b)) limits 
us to ramping times below a millisecond, such that we 
operate in the non-adiabatic regime, where significant 
Landau-Zener oscillations arc to be expected. To ac- 
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cess the adiabatic regime, we had to broaden the band 
gap to 0.8 E Iec by introducing larger deviations from C4- 
symmetry. In this case, we can maintain most atoms 
in the 4th band over the entire anti-crossing. As the 
anti-crossing is passed (with decreasing AV, i.e, from 
right to left in Fig. 1(d)) the condensate wave function 
undergoes a dramatic change. While on the right side 
of the anti-crossing a C4-invariant momentum spectrum 
with a leading zero momentum peak is observed, similar 
to that shown in Fig. 2(a), on the left side of the anti- 
crossing the spectrum acquires finite momentum charac- 
ter and reduced C2 symmetry with the leading Bragg- 
orders arising at ±(1, — l)hk, as shown in Fig. 4(a) for 
AV — AV C )/E ICC = —1.3. The rotational symmetries 
seen in the momentum spectra reflect those of the B4- 
band calculated for the respective values of AV. Corre- 
spondingly, on the right side of the anti-crossing the cal- 
culated zero momentum Bloch-function displays approx- 
imate C4-symmetry with local 3s-orbits and local ls- 
orbits in the deep and shallow wells, respectively, while 
on the left side, the deep wells comprise local 2p a! _ J/ -orbits 
aligned along the (1, — Indirection. 

The non-adiabatic case was studied for AE-i^a ~ 
0.26 -E rcc corresponding to 25 nK. After preparation of 
the condensate at AV = AVi, within 400 fis we tune to 
values AV < AV C on the left side of the avoided cross- 
ing, and record the time evolution of the momentum 
spectrum. An example for (AV^ — AV C )/E ICC = —3.0 
is shown in Fig. 4(b), where the visibility [28| of the 
±(1, — l)fiic-Bragg peaks (blue disks and green triangles) 
and that of the (0, 0)fifc-peak (red squares) are plotted 
versus the holding time after the jump over the anti- 
crossing. We observe significant oscillations, which can 
be qualitatively modeled as a beat between the single- 
particle Bloch-functions of the 4th and 2nd band at the 
T-point. A closer inspection shows, that the frequency 
spectra of these oscillations comprise small second har- 
monic components as is illustrated for the ( — l,l)hk- 
peak in Fig. 4(c), which is expected as a result of the 
non-linearity introduced by collisional interactions. The 
second harmonic contributions give rise to slightly in- 
creased (decreased) curvatures in the minima (maxima) 
of the oscillations in Fig. 4(b). Collisional relaxation is 
also responsible for the observed decay of the visibility. 
If an additional lattice potential along the z-direction 
is applied in order to increase the collision energy per 
particle, we find correspondingly decreased decay times. 
Evaluation of curves similar as in Fig. 4(b) for variable 
AV yields the blue disks in Fig. 4(d). The solid trace 
repeats the energy difference between the 4th and the 
2nd band from Fig. 1(d) using the measured value of e 
specified below Eq.(l). Despite neglecting collisions, the 
calculations without use of fitted parameters well approx- 
imate the observations. The small deviations of the ob- 
served frequencies towards values slightly below the sin- 
gle particle calculations cannot be reduced by choosing 



different values for rj, e x , e y in Eq.(l), but rather indicate 
the effect of collisions. A simplified non-linear Landau- 
Zener model with the three zero momentum Bloch func- 
tions of the band crossing in Fig. 1(c) as basis modes, a 
3x3 coupling matrix adjusted to reproduce the single- 
particle anti-crossing in Fig. 1(d), and a collision matrix 
accounting for the collision processes among the three ba- 
sis modes confirms the observed trend, however, without 
yielding quantitative agreement. 
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